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FOREWORD 


Contract  AF  04(6ll)-7^^8  between  the  Air  Force  Flight  Test  Center,  Edwards  Air 
Force  Base,  California,  Air  Force  Systems  Command,  United  States  Air  Force  and  the 
Research  Laboratories  of  United  Aircraft  Corporation  pertains  to  investigations  of 
the  characteristics  of  a  unique  gaseous-core  nuclear  rocket  concept.  These  investi¬ 
gations,  which  were  directed  toward  the  development  of  information  pertinent  to  the 
evaluation  of  the  feasibility  of  this  concept,  include  the  following  (l)  experimental 
determination  of  the  characteristics  of  vortex  flow,  (2)  design  and  fabrication  of 
equipment  to  provide  a  radial  temperature  gradient  in  a  vortex  flow,  (3)  theoretical 
calculations  of  the  characteristics  of  the  end-wall  boundary  layer  in  a  vortex  tube, 
and  (4)  theoretical  calculations  of  radiant  heat  transfer  in  the  proposed  engine. 

The  work  accomplished  under  the  contract  is  summarized  in  the  following  Research 
Laboratories'  reports  which  comprise  the  required  Final  Report: 

I.  Summary  Report  (Report  R-2k^k-k) 

II.  E.:perimental  Investigation  of  Characteristics  of  Confined  Jet-Driven 
Vortex  Flows  (Report  R-2494— 2) 

III.  Heat  Transfer  to  Confined  Vortex  Flows  by  Means  of  a  Radio-Frequency  Gas 
Discharge  (Report  R-2494— 3) 

IV.  Theoretical  Solutions  for  the  Secondary  Flow  on  the  End  Wall  of  a  Vortex 
Ttibe  (Report  R-2494-1,  Present  Report) 

V.  Methods  for  Calculating  Radiant  Heat  Transfer  in  High-Temperature  Hydrogen 
Gas  (Report  M-l492-l) 
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Report  R-249A-I 


Theoretical  Solutions  for  the  Secondeiry  Flow 
on  the  End  Wall  of  a  Vortex  Tube 
Contract  AF  04(6ll)-74L8 


SUMMARY 


Theoretical  solutions  are  presented  for  the  growth  of  the  tangential  and 
radial  components  of  both  laminar  and  turbulent  boundary  layers  In  th'  secondary 
flow  on  the  end  wall  of  a  vortex  flow  confined  by  a  cylindrical  tube  having  exit 
nozzles  located  on  the  axis  of  rotation  at  each  of  the  two  end  walls.  Solutions 
for  a  strong-vortex  case  (tangential  velocity  Inversely  proportional  to  radius) 
Indicate  that  the  weight  flow  entering  the  end-wall  boundary  layer  Is  proportional 
to  the  friction  between  the  tangential  flow  and  the  end  wall.  This  weight  flow 
Is  in  many  cases  greater  than  the  radial  weight  flow  outside  the  end-wall  boundary 
layer  required  to  maintain  a  strong  laminar  vortex  flow  (greater  by  a  factor  of 
approximately  3OO  for  a  turbulent  boundary  layer,  a  tube  length  of  10  times  the 
tube  radius,  and  a  Reynolds  number  based  on  tube  radius  and  on  the  tangential 
velocity  within  the  tube  of  10^).  Solutions  for  the  weak-vortex  case  in  which  the 
tangential  velocities  at  small  radii  are  less  than  that  for  the  strong-vortex  case 
are  obtained  by  allowing  for  the  reduction  in  vortex  strength  outside  the  end-wall 
boundary  layer  due  to  the  loss  of  radial  flow  to  the  end -wall  boundary  layer. 
Solutions  obtained  for  the  weak-vortex  case  indicate  that  some  of  the  weight  flow 
entering  the  end-wall  bo vr.dary  layer  near  the  outside  radius  of  the  tube  will 
return  to  the  main  flow  near  the  inside  radius  of  the  tube.  The  magnitude  of  the 
flow  within  the  end-wall  boundary  layer  offers  a  possible  explanation  of  the  weak- 
vortex  strength  which  has  been  measured  by  several  investigators  and  which  had 
heretofore  been  attributed  to  turbulence  within  the  primary  flow  outside  the  end- 
wall  boundary  layer. 

All  theoretical  solutions  obtained  are  based  on  assumed  shapes  for  the  tan¬ 
gential  and  radial  boiindary  layer  velocity  profiles  and  assumed  laws  relating  the 
tangential  and  radial  wall  friction  to  the  Reynolds  number  of  the  flow.  Sample 
calculations  for  the  variation  with  radius  of  tangential  velocity,  boundary  layer 
weight  flow,  and  boundary  layer  thickness  for  a  laminar  boundary  layer  profile 
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and  friction  law  are  given  for  both  the  strong-vortex  and  the  weak-vortex  cases, 
whereas  numerical  examples  for  a  turbulent  boundary  layer  profile  and  friction  law 
are  given  only  for  the  strong-vortex  case.  The  accuracy  of  results  calculated  from 
the  theory  should  be  limited  only  by  the  accuracy  of  the  velocity  profiles  and 
sheeir  formiilas  used  in  the  calculation  procedure. 

The  basic  theoretical  procedure  and  the  SEunple  calculations  for  laminar  boundary 
layers  were  developed  as  part  of  the  Research  Laboratories'  own  research  program. 

The  adaptation  of  the  theory  to  turbulent  boundary  layers  and  the  calculation  of 
the  examples  for  turbulent  boundary  layers  were  made  under  Contract  AF  04(6ll)-74L8 
with  the  Air  Force  Flight  Test  Center,  Edwards  Air  Force  Base,  California,  Air  Force 
Systems  Command,  U.  S.  Air  Force. 


CONCLUSIONS 


1 .  The  radial  mass  flow  in  the  boundary  layers  on  the  end  walls  of  a  cylindrical 
tube  containing  a  vortex  flow  is  proportional  to  the  friction  between  the  rotating 
flow  and  the  end  wall. 

2.  The  passage  of  flow  into  the  end-wall  boundary  layer  is  sufficient  in  many 
cases  to  cause  a  substantial  reduction  in  the  tangential  velocities  near  the  center- 
line  of  the  vortex  tube  from  those  which  would  be  obtained  without  an  end-wall 
boundary  layer.  This  calculated  reduction  in  vortex  strength  due  to  secondary  flow 
may  explain  the  loss  in  vortex  strength  which  has  been  measured  by  several  investi¬ 
gators  and  which  had  heretofore  been  attributed  to  turbulence  within  the  primary 
flow  outside  the  end-wall  boundary  layer. 

3.  If  the  magnitude  of  the  flow  into  the  boundary  layer  is  sufficient  to 
cause  a  reduction  of  vortex  strength,  some  of  this  boundary  layer  flow  will  return 
to  the  primary  flow  near  the  centerline  of  the  vortex  tube.  In  some  cases  the 
return  of  secondary  flow  to  the  primary  flow  region  is  sufficient  to  cause  a  radial 
outflow  in  the  primary  flow  at  certain  radii  within  the  vortex  tube. 
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INTRODUCTION 


Recently  there  has  been  a  renewed  interest  in  the  use  of  vortex  flows  for 
applications  in  magnetohydrodynamic  power  generation  (Ref.  1)  and  nuclear  rocket 
propulsion  (Ref.  2).  The  results  of  experimental  studies  of  confined  vortex 
flows,  such  as  those  described  in  Refs.  3>  and  indicate  that  these  flows 
are  considerably  more  complex  than  those  generally  assumed  in  one -dimensional 
analyses.  It  was  suggested  that  vortex  flows  are  tiurbulent  and  that  certain 
three-dimensional  phenomena,  which  are  not  well  understood,  might  have  a  primary 
effect  on  such  flows.  Attempts  at  developing  a  satisfactory  method  for  predict¬ 
ing  the  flow  patterns  obtained  with  a  confined  vortex  were  unsuccessful. 

A  further  understanding  of  the  characteristics  of  three-dimensional  vortex 
flows  was  reported  in  Ref.  6.  In  particular,  the  existence  of  a  two-cell  vortex 
in  an  open-ended  vortex  tube  was  demonstrated.  However,  these  results  are  not 
applicable  to  vortex  flows  confined  by  end  walls.  Another  analysis,  reported  in 
Ref.  7,  considered  vortex  flows  in  cylindrical  tubes  having  end  walls  and  included 
a  friction  factor  to  account  for  the  end-wall  boundary  layer.  However,  because 
this  analysis  was  also  one-dimensional,  the  results  cannot  be  used  to  predict  any 
of  the  three-dimensional  effects  existing  in  vortex  tubes . 

Finally,  subsequent  to  the  completion  of  current  studies  on  vortex  flows  at 
the  UAC  Research  Laboratories,  two  papers  on  such  flows  (Refs.  8  and  9)  were 
received  which  apply  directly  to  the  work  reported  herein.  Reference  8  is  an 
analytical  study  of  the  laminar  boundary  layer  on  a  disc  in  a  rotating  flow  with 
prescribed  tangential  velocity  distributions.  The  analytical  procedure  is  essen¬ 
tially  that  used  in  this  report.  However,  because  this  analysis  considered  only 
the  semi-infinite  vortex,  the  complex  interaction  between  the  primary  and  secondary 
flows  was  not  examined.  Reference  9  is  an  experimental  study  of  vortex  flow 
which  indicates  that  as  much  as  905t  of  the  vortex  mass  flow  may  be  carried  in  the 
end -wall  boundary  layer. 

The  object  of  the  study  reported  herein  was  to  formulate  a  three-dimensional 
theory  for  confined  vortex  flows  in  order  to  calculate  the  effects  of  the  end-wal 
boundary  layer  and  secondary  axial  flow  on  the  distribution  of  angular  momentum 
and  tangentisd.  velocity  in  the  primary  flow  of  the  vortex. 

\ 
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FLOW  MODEL  AND  BASIC  ASSUMPTIONS 


A  diagram  of  the  analytical  model  used  in  the  analysis  is  shown  in  Fig.  1. 

The  vortex  is  assumed  to  be  contained  in  a  cylinder  of  radius  ft  and  length  z  i  , 
and  the  flow  is  uniformly  injected  into  the  cylinder  at  ri  with  a  tangential 
velocity  ,  and  a  radial  velocity  urt  •  Two  exit  nozzles  of  radius  r,  are 
located  at  each  end  of  the  vortex  tube  concentric  with  the  axis  of  the  cylinder. 

A  sketch  showing  the  velocities  considered  in  the  analysis  is  given  in  Fig.  2. 
The  secondary  flow  in  the  end-wall  boundary  layer  is  characterized  by  large  radial 
velocities  relative  to  the  radial  velocities  in  the  primary  flow.  The  origin 
of  this  radial  flow  can  be  understood  from  the  following  discussion.  In  the  pri¬ 
mary  flow  the  centrifugal  forces  produced  by  the  tangential  motion  of  the  fluid 
are  balanced  by  the  radial  pressure  gradient.  However,  in  the  end-wall  boundary 
layer  the  tangential  velocity  is  reduced  to  zero  at  the  wall  by  shear  forces,  but 
the  radial  pressure  gradient  remains.  This  radial  pressure  gradient  accelerates 
the  flow  in  the  boundary  layer  to  large  radial  velocities  relative  to  the  radleQ. 
velocities  in  the  primary  flow.  Thus,  although  the  boundary  layer  thickness  is 
small,  an  appreciable  fraction  of  the  radial  flow  may  pass  through  this  end-wall 
boundary  layer. 

The  following  assumptions  are  employed  in  the  analyses  discussed  in  this 
report ; 

1.  The  flow  is  incompressible. 

2.  The  flow  is  rotationally  symmetric. 

3.  The  pressure  gradient  in  the  axial  direction  in  both  the  primary 
and  secondary  flows  is  zero. 

4.  The  thickness  of  the  boundary  layer  on  the  end  walls  is  small 
relative  to  the  radius  of  the  vortex  tube. 

5.  The  thickness  of  the  boundary  layer  in  the  tangentieO.  direction  is 
equal  to  the  thickness  of  the  boundary  layer  in  the  radial  direction. 

6.  Both  the  radial  and  tangential  velocity  profile  shapes  within  the 
boundary  layer  are  independent  of  radius. 

7.  The  friction  governing  the  shear  between  the  flow  within  the  boundary 
layer  and  the  end  wall  of  the  tube  is  a  function  of  the  Reynolds  number 
based  on  the  thickness  of  the  boundary  layer  in  the  end  wall.  (See 
Appendix  A ) . 
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8.  The  flow  outside  the  boundary  layer  is  laminar  (see  discussion 

in  section  entitled  Limit  of  Applicability  of  Strong -Vortex  Solut  on). 

9.  The  radial  velocity  is  much  less  than  the  tangential  velocity  in  the 
primary  flow  outside  the  end -wall  boundary  layer. 

10.  The  tangential  velocity  outside  the  end-wall  boundary  layer  is  inde¬ 
pendent  of  axial  distance.  This  assumption  is  Justified  by  the  data 
presented  in  Refs.  1  and  10 and  is  proven  analytically  in  Appendix  A 
for  the  special  case  with  the  tangential  velocity  much  larger  than 
the  radial  velocity. 

11.  The  radial  velocity  outside  the  end-wall  boundary  layer  is  independent 
of  axial  distance. 

12.  The  axial  velocity  of  the  flow  passing  through  the  exit  nozzles  (see 
Fig.  1)  is  independent  of  radius  within  these  nozzles. 


GENERAL  FORM  OF  SOLUTIONS 


I 

1. 

I* 


1: 


The  general  form  of  the  solutions  to  be  obtained  can  be  shown  by  considering 
a  torque  balance  at  the  end  wall  of  a  vortex  tube.  In  this  torque  balance  the 
moment  created  by  the  shear  force  between  the  moving  fluid  and  the  wall  must  be 
coimteracted  by  a  moment  created  by  the  change  in  tangential  velocity  of  the  flow 
entering  the  boundary  layer. 


where 


Ws 


is  mass  flow  entering  boundary  layer. 

is  tsingential  velocity  -  proportional  to  change  in  tangential 
velocity  of  flow  entering  boundary  layer. 


(1) 


ri 

is 

Cf 

is 

P 

7rr, 

is 

is 

^^<^1  is  dynamic  pressure  of  rotating  gas . 
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Therefore,  the  mass  flow  within  the  boundary  layer  on  one  end  wall  of  a  vortex 
tube  is  proportional  to  the  right-hand  side  of  the  following  equation. 

However,  from  the  definition  of  tangential  Reynolds  number 


Re 


f 


(3) 


Wg  ~  Re,  (4) 

It  is  assumed  that  the  friction  coefficient  varies  with  Reynolds  nusiber  the 
same  as  on  a  flat  plate .  For  laminar  flow 


I 


Therefore, 


SL 


K 

L 


{2-nT|/i.)  Re 


1/2 

f 


(5) 


(6) 


The  coefficient  Kl  in  Eq.  (6)  must  be  chosen  to  account  for  all  of  the  constants 
missing  in  the  preceding  equations. 


For  turbulent  flow 


I 


(7) 


Therefore, 


W. 


ST 


K^(27rr,  ^)Re 
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It  is  also  assvmed  that  the  eaqiressions  for  the  veuriation  of  boundary  layer 
thickness  with  Reynolds  number  are  similar  to  the  corresponding  expressions  for 
flow  over  a  flat  plate . 


(9) 


Ct 


Re 


1/5 


(10) 


The  procediires  employed  to  develop  relations  for  the  four  coefficients,  Kl  , 

Kt  ,  ,  and  which  are  required  to  evaluate  Eqs.  (6),  (8),  (9),  and  (lO) 

are  developed  in  Appendixes  A  and  B. 


STRONG-VORTEX  SOLOTION 


The  results  discussed  in  the  following  section  were  obtained  on  the  basis 
that  the  vortex  is  sufficiently  strong  that  the  angular  momentum  of  the  flow 
outside  the  boundary  layer  is  constant  (velocity  outside  the  boundary  layer 
inversely  proiwrtional  to  radius). 


Specific  Solutions  for  Strong-Vortex  Case 

A  coiiq>lete  solution  of  the  flow  in  the  boundary  layer  on  the  end  wall  of  the 
vortex  tube  requires  consideration  of  the  radial  momentvun  equation  as  well  as  the 
tangential  momentum  equation  employed  in  the  preceding  section.  In  this  radlsG. 
momentum  eqiiation,  the  radial  acceleration  of  the  flow  within  the  boundary  layer 
is  determined  by  the  forces  resulting  from  the  radial  static  pressure  gradient, 
the  radial  shear  on  the  end  wall  of  the  vortex  tube,  and  the  change  in  radial 
velocity  of  the  mass  flow  entering  the  boundary  layer  from  the  primary  flow  outside 
the  boundary  layer.  Procedures  for  simultaneous  solution  of  the  tangential  and 
radial  momentum  equations  are  discussed  in  Appendix  A. 

To  obtain  a  solution  for  the  tangential  euid  radial  momentum  equations,  it  is 
necessary  to  assume  a  shape  for  both  the  tangential  and  radial  velocity  profiles. 
Plots  showing  the  tangential  and  radial  velocity  profiles  which  have  been  employed 
in  the  analyses  are  shown  in  Figs.  3  and  A.  Two  of  the  three  profiles  shown  in 
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each  of  the  figures  eu:e  for  laminar  flow,  whereas  the  third  profile  shape  is  for 
turbulent  flow.  Results  of  calculating  the  value  of  Kl  at  the  exit  radius  for 
different  combinations  of  laminar  velocity  profiles  are  shown  in  Table  I.  It  can 
be  seen  from  Table  I  that  the  value  of  Kl»  (which  is  employed  in  calculating  the 
weight  flow  in  the  boundary  layer  in  Eq.  (6)  at  R  =  R*  )  varies  by  a  factor  of 
slightly  more  than  two  for  the  two  extreme  results .  In  an  effort  to  choose  the 
most  likely  laminar  boundary  layer  profile  to  be  encountered  in  a  vortex  tube, 
each  of  the  profile  combinations  shown  in  Table  I  was  used  in  an  analysis  of  the 
secondary  flow  on  a  rotating  disc  (see  Appendix  B).  The  results  of  con^arlng  the 
boundary  layer  mass  flow  calculations  on  a  rotating  disc  with  a  known  exact  solu¬ 
tion  indicated  that  velocity  profiles  and  (^  in  Figs.  3  4  would  result 

in  the  best  agreement  with  the  exact  solution.  (The  secondary  mass  flow  on  a 
rotating  disc  depends  directly  on  the  parameter  B.  Table  II  presents  a  comparison 
of  the  exact  value  of  B  with  the  approximate  values  given  by  the  several  assumed 
velocity  profiles.)  According  to  Table  I,  the  value  of  Kl*  for  profiles  and 

@ 


Xl. 


1.226 


(11) 


Therefore  the  mass  flow  in  the  boundary  layer  on  one  end  wall  of  a  vortex  tube 
at  a  v£j.ue  of  R*  of  0.1  is,  from  Eq.  (6), 

1/2 

'^SL*  ■  '  226  (2  irr,  At)  Re,  (12) 

The  value  of  Kt*  for  turbulent  flow  for  profiles  and  given  in  Figs. 

3  and  4  is 


K_  =  O.H3 

T  • 


(13) 


The  corresponding  expression  for  the  mass  flow  in  the  boundary  layer  on  one  end 
wall  of  a  vortex  tube  at  a  value  of  R,  of  0.1  is,  from  Eq.  (8), 


W. 


ST* 


0.113  (2Trr,  At)Re 


4/6 

t 


(1^) 
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Curves  determined  from  Eqs.  (12)  and  (lA)  are  given  as  a  function  of  tangen¬ 
tial  Reynolds  number  in  Fig.  5*  As  noted  on  the  figure^  the  transition  Reynolds 
number  observed  in  tests  of  rotating  discs  is  on  the  order  of  10^  (see  Ref.  9)* 

A  plot  showing  the  variation  with  radius  of  the  dimensionless  radial  mass 
flow  within  the  boundary  layer  on  one  end  wall  of  a  vortex  tube  is  given  in 
Fig.  6  for  both  the  laminar  flow  and  turbulent  flow  solutions.  The  values  given 
in  Fig.  6  for  a  value  of  R,  =  0.1  correspond  to  the  constants  in  Eqs.  (ll) 
through  (lA). 

An  analytical  solution  is  given  in  Appendix  A  for  the  variation  of  boundary 
layer  thickness  with  radius  for  a  laminar  boundary  layer  on  the  end  wall  of  a 
vortex  tube.  The  dimensionless  boundary  layer  thickness  for  R  between  R*  euid 
1.0  from  Eq.  (l22)  of  Appendix  A  is  as  follows: 


,  ,1/4  2/S  ,  4/5, 1/4 

0.93l(-0|)  R  (l-R  ) 


The  constant  in  the  preceding  equation  is  presented  in  Table  I  for  four 

different  combinations  of  laminar  velocity  profiles .  It  can  be  seen  from  Table  I 
that  the  veuriation  of  this  constant  among  the  different  profiles  is  less  than 
the  variation  of  Kl,  among  the  seune  profiles.  The  variation  of  Cl  with  radius 
for  laminar  profiles  and  in  Figs.  3  suid  4  is  given  in  Fig.  7-  The  varia¬ 

tion  of  Cl  with  radius  for  the  different  laminar  profiles  in  Tuble  I  is  identical. 
However,  the  absolute  value  of  Cl*  any  value  of  R,  is  proportional  to  the 
value  of  the  constant  in  Table  I. 


Also  shown  in  Fig.  7  Is  the  variation  of  Ct  with  radius  for  a  turbulent 
boundary  layer  having  profile  shapes  and  of  Figs.  3  and  4.  The 

numerical  solution  used  to  obtain  Ct  Tor  the  turbulent  layer  is  described  in 
Appendix  A. 


Limit  of  Applicability  of  Strong-Vortex  Solution 

The  strong-vortex  solutions  are  valid  as  long  as  the.e  is  sufficient  radial 
mass  flow  in  the  primary  flow  outside  the  end-wall  boundary  layer  to  provide  the 
tangential  momentum  to  overcome  the  shear  forces  which  exist  near  the  centerline 
of  the  vortex.  This  required  radial  mass  flow  is  usually  given  in  terms  of  the 
radial  Reynolds  number,  RCf ,  defined  as  follows: 


pu  r,  W 

I  _  _ 

p.  '  ZTT/XZ, 


(16) 
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The  effect  of  radial  Reynolds  number  on  the  variation  of  tangential  velocity 
with  radius  for  the  case  of  zero  secondary  flow  is  given  in  Fig.  8.  It  can  be 
seen  from  Fig.  8  that  the  vortex  strength  begins  to  fall  off  rapidly  between  a 
radial  Reynolds  number  of  10  and  5*  The  criterion  given  in  Fig.  8  does  not  hold 
when  there  is  a  substantial  amount  of  flow  within  the  end-wall  boundary  layer. 

For  this  case,  an  effective  Reynolds  number  based  on  the  flow  outside  the  end- 
wall  boundary  layer  must  be  employed.  If  a  is  defined  as  the  ratio  of  the  flow 
within  the  boundary  layer  on  both  end  walls  of  the  vortex  tube  to  the  total  flow 
passing  through  the  vortex  tube,  then. 


(17) 


In  the  preceding  equation  O’*  is  the  value  of  O’  at  the  exit  radius  where  the 
greatest  amount  of  flow  has  entered  the  boundary  layer.  From  the  definition 
of  O’  , 


O’ 


2W, 


2irr^Z|  pUf 


(18) 


Substituting  from  Eq.  (6),  which  gives  the  flow  in  the  end-wall  boundary  layer 
with  a  laminar  velocity  profile. 


2K  (TTf 

O’  =  - 4 - ! - L_ 

27rr,z, 


(19) 


From  Eq.  (l6) 


Since 


O’l  =  2Kl 


Re 


i/z 


Re, 


hi 

Re^ 


(20) 


(21) 
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Therefore, 


2K^  ^ 


-L  =  ^  -^)  =  \-8. 

\  RCj  Uj.  z,  / 

where  the  secondary  flow  parameter  for  laminar  flow,  ;9|_ 

2 


= 


^  Ur, 


/3u  =  2 


Re 


1/2 


Re  r 


In  the  same  manner,  for  turbulent  flow. 


/  2r, 


Ret 


4/5 


Re , 


where  the  secondary  flow  parameter  for  turbulent  flow,  /; 

2 


= 


Re 


"Vs  0 


t 
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(23) 

,  is  defined  as  follows: 

(24) 

(25) 

(26) 

(27) 

It  ,  Is  given  below. 

(28) 
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R.y’ 


Substituting  Eqs.  (23)  and  (27)  into  Eq.  (l?)  yields 


(^^r)  «ff  ‘ 
Re^ 


( r )  eff 

Rer  . 


(29) 


(30) 


(31) 


As  noted  in  a  preceding  paragraph^  an  effective  radial  Reynolds  nunber  of 
between  5  and  10  is  required  to  provide  a  strong  vortex  outside  the  end-wall 
boundary  layers.  Therefore,  values  of  (Rcr)  ^ff  of  5  and  10  were  used  in 
evaluating  Eqs.  (30)  and  (3l)»  The  results  of  this  evaluation  are  presented  in 
Fig.  9*  It  can  be  seen  from  Fig.  9  that  the  maximum  permissible  value  of  the 
secondary  flow  parameter  for  applicability  of  the  strong-vortex  results  increases 
with  an  increase  in  radial  Reynolds  number. 

The  limit  of  applicability  for  the  strong-vortex  results  can  also  be  presented 
in  terms  of  tangential  Reynolds  number  by  combining  Eqs.  (30)  and  (31)  with 
Eqs.  (25)  and  (29)-  The  results  obtained  in  this  manner  using  a  value  of(Rer)eff 
of  5  are  presented  in  Fig.  10  for  three  different  values  of  Z|/r,  .  The  tangen¬ 
tial  Reynolds  number  must  be  lower  than  that  shown  in  Fig.  10  if  the  strong- 
vortex  results  are  to  apply.  The  results  presented  in  Figs.  9  and  10  were  derived 
on  the  basis  that  the  vortex  strength  was  determined  entirely  by  the  radied.  mass 
flow  at  a  radius  equal  to  .  In  practice,  the  tangential  velocity  profile  is 
determined  by  some  kind  of  average  of  the  radial  mass  flows  which  exist  oui^side 
the  bour.dary  layer  at  all  values  of  radius .  A  more  exact  calculation  of  the  limit¬ 
ing  value  of  /3i_  for  applicability  of  the  laminar  strong-vortex  results  will  be 
described  in  a  following  section. 

Data  on  vortex  strength  are  presented  in  Refs.  3^  and  5  for  radial  Reynolds 
numbers  ranging  from  20  to  2000.  According  to  Fig.  8  these  radial  Reynolds  numbers 
should  produce  strong  vortices  if  the  fraction  of  the  flow  in  the  end-wall  boxindary 
layers  is  small  and  if  the  flow  outside  the  end-wall  bovuidary  layers  is  laminar. 
However,  the  velocity  profiles  reported  in  Refs.  3  and  4  corresponded  to  radial 
Reynolds  numbers  on  the  order  of  2  to  4.  These  differences  between  measiored 


CONFIDENTIAL 


12 


CONFIDENTIAL 

Reynolds  numbers  based  on  total  mass  flow  and  Reynolds  numbers  determined  from 
velocity  profiles  were  attributed  in  Refs.  3  and  4  to  turbulence  outside  the 
end-wall  boundary  layer.  However,  this  same  result  could  be  obtained  by  the 
passage  of  a  large  percentage  of  the  total  radial  mass  flow  through  the  end-wall 
boundary  layers.  To  check  this  possibility,  the  results  presented  in  Fig.  5 
were  employed  to  calculate  the  secondary  flow  for  the  data  presented  in  Refs. 

3,  4,  and  5*  The  resulting  calculated  values  of  secondary  flow  were  divided 
by  the  total  measured  mass  flow  in  each  test  and  the  results  are  presented  in 
Fig.  11  as  a  function  of  tangential  Reynolds  number.  It  can  be  seen  from 
Fig.  11  that  the  calculated  secondary  flows  are  higher  than  the  total  mass  flow 
in  the  vortex  tubes  for  most  of  the  test  points .  Therefore,  it  is  possible  that 
all  of  the  loss  in  vortex  strength  reported  in  Refs.  3  and  4  is  due  to  secondary 
flow  rather  than  to  turbulence .  The  significance  of  values  of  secondary  flow 
calculated  by  the  strong-vortex  theory  being  greater  than  the  total  flow  passing 
through  the  vortex  tube  will  be  discussed  in  the  following  section.  The 
hypothesis  of  loss  in  vortex  strength  due  to  secondary  flow  is  also  strengthened 
by  a  measurement  reported  in  Ref.  j  in  which  90?^'  of  the  total  radial  mass  flow 
in  a  vortex  tube  was  found  to  pass  down  the  end-wall  boundary  layer.  Insuffi¬ 
cient  data  were  reported  in  Ref.  9  to  permit  calculation  of  the  theoretical 
weight  flow  corresponding  to  this  measured  weight  flow. 


WEAK-VORTEX  SOLLfTION 


As  noted  in  the  preceding  section,  the  strong-vortex  solution  is  no  longer 
valid  when  the  mass  flow  passing  into  the  end-wall  boundary  layer  is  so  great 
that  there  is  insufficient  radial  flow  outside  the  boundary  layer  to  maintain  a 
strong  vortex.  The  solution  of  the  weak-vortex  case  requires  simultaneous  solu¬ 
tion  of  the  boundary  layer  and  primary  flow  equations .  The  procedure  employed 
is  described  in  Appendix  A.  All  calculations  discussed  in  the  following  sub¬ 
sections  were  obtained  on  the  basis  of  laminar  flow  in  the  boundary  layers  with 
profiles  and  in  Figs.  3  and  4.  However,  the  qualitative  nature  of  the 

results  is  expected  to  apply  for  turbulent  boundary  layers . 


Solutions  for  Fixed  Radial  Reynolds  Number 

A  series  of  solutions  was  obtained  for  a  fixed  value  of  the  radial  Reynolds 
number  of  10  and  for  a  series  of  values  of  the  secondary  flow  parameter,  /Sl  . 

It  can  be  seen  from  Eq.  (25)  that  an  increase  in  (i\_  for  a  fixed  vortex  tube 
geometry  and  a  fixed  value  of  radial  Reynolds  number  represents  an  increase  in 
tangential  Reynolds  number.  The  effect  of  an  increase  in  /3l  on  the  variation 
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of  angular  momentiim  and  tangential  velocity  with  radius  is  shown  on  Figs.  12 
and  13.  It  can  be  seen  from  these  figures  that  an  increase  in  and  the 

accompanying  increase  in  secondary  flow  (see  Eqs.  (l?)  and  (23)  )  results  in  a 
decrease  in  vortex  strength  near  the  centerline  of  the  vortex  tube.  The  varia¬ 
tion  of  static  pressure  coefficient  with  radius  is  shown  in  Fig.  l4  for  several 
values  of  (i\_  .  Since  the  change  of  pressure  with  radius  is  directly  proportional 
to  the  square  of  the  tangential  velocity,  the  differences  in  static  pressure 
distributions  shown  in  Fig.  l4  for  the  different  values  of  (i\_  £u:e  directly 
due  to  the  decrease  of  vortex  strength  with  increasing  /?l •  The  curves  in 
Figs.  13  and  l4  for  =  0  are  identical  to  those  employed  in  the  analyses  in 
Refs.  3j  4,  and  5- 

The  radial  mass  flow  in  the  end-wall  boundary  layer  is  shown  as  a  function 
of  radius  in  Figs .  15  and  16  for  different  values  of  the  secondary  flow  parameter, 
/?!_  .  It.  can  be  seen  from  Fig.  15  that  the  nondimens ional  mass  flow  within  the 
end-wall  boundary  layer  at  small  values  of  radius  decreases  with  an  increase 
in  /3l  ;  i.e.,  with  a  decrease  in  vortex  strength  near  the  centerline  of  the  tube 
(see  Figs.  12  and  13)^  although  the  fraction  of  secondary  flow  increases.  It  can 
also  be  seen  from  Figs .  15  and  I6  that  the  radial  mass  flow  in  the  end-wall  boun¬ 
dary  layer  decreases  with  decreasing  radius  near  the  centerline  of  the  vortex  tube. 
The  reason  for  this  reversal  of  trend  in  the  variation  of  boundary  layer  mass  flow 
with  radius  is  as  follows.  Near  the  outer  radius  of  the  tube  there  is  more  of  a 
tendency  to  maintain  constant  angular  momentijm  in  the  primary  flow  than  in  the 
secondary  flow;  therefore,  flow  passes  from  the  primary  to  the  secondary  in  an 
effort  to  maintain  angular  momentum  in  the  secondary  flow.  However,  at  low 
values  of  radius  the  tendency  to  maintain  angular  momentum  in  the  secondary  flow 
is  greater  than  in  the  primary  flow;  therefore,  there  is  a  tendency  for  flow  to 
pass  from  the  secondary  flow  to  the  primary  flow  in  an  effort  to  maintain  a  strong 
vortex  in  the  primary  flow.  The  significance  of  values  of  cr  greater  than  unity 
in  Fig.  16  is  discussed  in  a  following  paragraph.  The  passage  of  flow  from  the 
boundary  layer  to  the  main  stresim  is  also  predicted  by  the  theory  of  Ref.  8  and 
some  experimental  verification  is  given  by  the  ink-dye  experiments  reported  in 
Ref.  9. 

The  variation  of  dimensionless  boundary  layer  thickness  with  radius  for 
several  values  of  the  secondary  flow  parameter,  (i\_  ,  is  shown  in  Fig.  17 .  It 
can  be  seen  from  Fig.  17  that  an  increase  in  results  in  an  increase  in 
dimensionless  boundary  layer  thickness  even  though  it  resiilts  in  a  decrease  in 
mass  flow  within  the  boundary  layer. 

The  streamlines  formed  by  axial  and  radial  velocities  for  a  particular 
vortex  tube  geometry  are  plotted  in  Figs .  18  throu^  20  for  values  of  the 
secondary  flow  parameter  of  0.4,  1.0  and  4.0  and  for  Re^  =  10.0.  The 
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geometry  of  the  vortex  tube  Is  defined  by  =  0.10  and  Z,  /r,  =  6.0.  Since 

the  vortex  tube  is  symmetrical  about  the  axis  and  the  z/z,  =  l/2  plane,  the 
streamlines  are  plotted  for  only  one  quadrant.  The  tangential  velocity  distri¬ 
bution  is  also  shown  on  these  figures.  Of  particular  interest  are  the  streamlines 
in  Fig.  20  which  illustrate  the  "vortex  cell"  condition  in  which  there  is  radial 
outflow  in  the  primary  flow.  In  order  to  preserve  continuity,  the  inward  flow 
in  the  end-wall  boundary  layer  is  equal  to  the  sum  of  the  outward  flow  in  the 
vortex  cell  and  the  total  mass  flow  passing  through  the  vortex  tube.  Ihis  vortex 
cell  condition  arises  when  the  secondary  flow  necessary  to  overcome  friction  on 
the  end  wall,  is  greater  than  the  flow  injected  into  the  vortex  tube. 


Effect  of  Radial  Reynolds  Number  on  Weak-Vortex  Solution 

The  variation  of  angular  momentum,  tangential  velocity,  secondary  mass  flow, 
and  boundary  layer  thickness  with  radius  for  several  values  of  the  secondary  flow 
parameter,  ,  are  presented  in  Figs.  21  through  24  for  a  radial  Reynolds  nvunber 
of  5  and  in  Figs.  25  through  28  for  a  radial  Reynolds  number  of  20.  It  can  be 
seen  by  comparing  Figs.  21  and  22  with  Figs.  12  and  13  that  a  decrease  in  radial 
Reynolds  number  results  in  an  increased  sensitivity  of  tangential  velocity  profiles 
to  the  secondary  flow  parameter,  /9l  -  The  converse  can  be  seen  by  conqparing  Figs. 
25  and  26  with  Figs.  12  and  13 . 

The  ratios  of  secondary  flow  to  total  flow  obtained  from  Figs.  l6,  23,  and 
27  are  summarized  in  Figs.  29  and  30*  The  maximum  value  of  c  ,  denoted  as  <r|^  , 
is  presented  in  Fig.  29  whereas  the  value  of  cr  at  the  nozzle  exit  radius,  cr^  , 
is  presented  in  Fig.  30.  It  ceui  be  seen  from  these  figures  that  the  secondary 
flow  ratio  decreases  from  that  calculated  for  the  strong-vortex  solution  as  the 
secondary  flow  parameter  is  increased.  The  fraction  of  flow  recirculated;  i.e., 
the  flow  which  is  returned  from  the  boundary  layer  to  the  primary  flow  divided  by 
the  total  flow  passing  through  the  vortex,  <7^-  cr^  ,  is  presented  in  Fig.  31. 

These  results  were  obtained  from  Figs.  29  and  30. 

The  preceding  discussions  indicate  the  existence  of  several  flow  regimes  in 
a  vortex  tube.  Tnese  flow  regimes  are  shown  in  Fig.  32.  Two  separate  criteria 
are  used  to  separate  the  different  flow  regimes:  one  criterion  describing  the 
strength  of  the  vortex  (solid  lines  in  Fig.  32),  and  the  other  criterion  describing 
the  fraction  of  the  flow  passing  through  the  end-wall  boundary  layer  (dashed  lines 
in  Fig.  32).  The  curve  separating  the  strong-vortex  regime  from  the  weak-vortex 
regime  is  defined  such  that  the  tangential  velocity  at  the  exit  radius  is  equal 
to  95?^  of  the  tangential  velocity  which  would  exist  in  an  inviscid  vortex.  (This 
curve  on  Fig.  32  falls  midway  between  the  two  curves  for  laminar  flow  given  on 
Fig.  9*)  The  curve  separating  the  weak-vortex  regime  from  the  very  weak-vortex 
regime  is  defined  such  that  the  tangential  velocity  at  the  exit  radius  is  equal 
to  the  tangential  velocity  at  the  outside  radius  of  the  vortex  tube. 
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The  secondary  flow  regimes  are  distinguished  by  the  values  of  cTj^  ,  the 
ratio  of  the  maximvim  value  of  the  flow  in  the  end -wall  boundary  layers  to  the 
total  flow  passing  through  the  vortex  tube.  The  small  secondary  flow  regime  is 
defined  such  that  is  less  than  0.1.  The  large  secondeiry  flow  regime  is 

defined  where  is  between  0.1  and  1.0.  Finally,  the  vortex  cell  regime  is 

defined  where  is  greater  than  1  such  that  some  of  the  fluid  is  trapped  in  a 

vortex  cell,  as  shown  by  the  streamline  pattern  in  Fig.  20.  It  can  be  seen  from 
Fig.  32  that  there  is  a  regime  where  there  are  large  secondary  flows  and  a  strong 
vortex  and  also  a  regime  where  there  are  trapped  vortex  cells  and  a  weak  vortex 
motion. 


CONFIDENTIAL 


R. 2494-1 


CONFIDENTIAL 

REFERENCES 


1.  Donaldson,  C.  duP.:  The  Magnetohydrodynamic  Vortex  Power  Generator,  Basic 
Principles  and  Practical  Problems.  ARAP  Report  No.  30^  1961* 

2.  Kerrebrock,  J.  L.  and  R.  V.  Megreblian:  An  Analysis  of  Vortex  Tubes  for 
Combined  Gas  Phase  Fission-Heating  and  Separation  of  the  Fissionable  Material. 
ORNL  CF- 57-11- 3  Rev.  1,  April  1958. 

3.  Keyes,  J.  J.:  An  Experimental  Study  of  Gas  Dynamics  in  High  Velocity  Vortex 
Flow.  Proc.  Heat  Transfer  and  Fluid  Mechanics  Inst.,  Stanford  University, 

June  i960. 

4.  Keyes,  J.  J.,  and  R.  E.  Dial:  An  Experimental  Study  of  Vortex  Flow  for 
Application  to  Gas-Phase  Fission  Heating.  0RNL-2837>  April  i960. 

5.  Ragsdale,  R.  G.:  NASA  Research  on  the  Hydrodynamics  of  the  Gaseous  Vortex 
Reactor.  NASA  TN  D-288,  September  i960. 

6.  Donaldson,  C.  duP.  and  R.  D.  Sullivan;  Behavior  of  Solutions  of  the  Navier- 
Stokes  Equations  for  a  Complete  Class  of  Three-Dimensional  Viscous  Vortices, 
Proc.  Heat  Trans,  and  Fluid  Mech.  Inst.,  Stanford  University,  June  i960. 

7.  Rietema,  K.,  and  H.  J.  Krajenbrink:  Theoretical  Derivation  of  Tangential 
Velocity  Profiles  in  a  Flat  Vortex  Chamber  -  Influence  of  Turbulence  and 
Wall  Friction.  App.  Sci.  Res.,  Sec.  A,  Vol.  8,  No.  223,  1959* 

8.  Mack,  L. :  Laminar  Boundary  Layer  on  a  Disk  in  a  Rotating  Flow.  JPL  Research 
Summary  No.  36-9^  July  1,  1961. 

9.  Kendal*,  J.:  Vortex  Separation  Experiments.  JPL  Research  Summary  No.  36-3/ 
June  15,  i960. 

10.  Hartnett,  J.  P.  and  E.  R.  G.  Eckert:  Experimental  Study  of  the  Velocity  and 
Temperature  Distribution  in  a  High  Velocity  Vortex  Tube.  Proc.  Heat  Transfer 
and  Fluid  Mech.  Inst.,  1956. 

11.  Schlichting,  H.:  Boundary  Layer  Theory.  Pergamon  Press,  1955- 

12.  Hall,  M.  S.:  The  Momentvun  Integral  Equations  for  Three-Dimensional  Laminar 
Boundary  Layers  in  Incompressible  Flow.  Australian  Aeronautical  Research 
Committee  Report  ACA-62.  November  1959* 


CONFIDENTIAL 


17 


18 


r.2494-1 


CONFIDENTIAL 

LIST  OF  SYMBOIS  (Cont.) 


Tangential  Reynolds  number. 


Reynolds  number  based  on  boundary  layer  thickness. 


p  0,^,8 


Torque  per  unit  radius,  lb 

Radial,  tangential,  axial  velocities,  ^t/sec 

Radial,  tangential,  axial  velocities  In  primary  flow,  ft/sec 

Ur 

Up  Dimensionless  radleil  velocity  In  primary  flow,  -= - 

u<^ 

Dimensionless  tangential  velocity  in  primary  flow,  q ^ 
g  Axial  velocity  entering  boundary  layer 

Wc  , 

Vp  Average  radial  velocity  in  boundary  layer,  2TTp^\j — 

Wg  Secondary  mass  flow  in  one  end-wall  boxmdary  layer,  slugs/sec 

z  Axial  distance,  ft 

Z  Axial  distance  in  bovindary  layer  (dimensionless),  z/S 

Secondary  flow  parameter  (laminar),  ,,3  Yj- 

2  r, 

Pt  Secondary  flow  parameter  (turbulent),  — — ^  7; —  -= — 

'  Re| u  r I  4 1 

8  Boundary  layer  thickness 

€  Small  quantity  (dimensionless) 

Dimensionless  boundary  layer  thickness  for  a  vortex  tube  (laminar  - 
S/r,  )(turbulent  -  Re  j  8/  r,  ) 

1/2 

17  Dimensionless  boundary  layer  thickness  for  rotating  disc  (  Re  t  S/r,  ) 

K  Dimensionless  singular  momentum,  — = - 

'  I <^i 

fj.  Molecular  viscosity,  lb  sec/ft^ 

V  Kinematic  viscosity,  ft^/sec 
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LIST  OF  SYMBOIS  (Cont.) 


P 


^4> 


CJ 


Dummy  variable  of  integration 
Density,  slugs/ft ^ 

2Ws 

Secondary  flow  ratio,  -tt— - = - 

'  27rr,pur^z, 

Radial  shear  stress,  Ib/ft^ 

Tangential  shear  stress,  Ib/ft^ 

Stream  function 

\1/ 

Dimensionless  stream  function,  - =— 

r,  z,Ur 

Angular  velocity  of  rotating  disc 


I 


Subscripts 

I  Quantity  evaluated  at  outer  radius  of  vortex  tube 

*  Quantity  evaluated  at  exit  radius  of  vortex  tube 

L  Quantity  evaluated  for  laminar  boundary  layers 

M  Maximum  value  of  quantity 

T  Quantity  evaluated  for  turbulent  boundary  layers 
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APPENDIX  A 


DERIVATION  AND  SOLUTION  OF  THE  EQUATIONS 
DESCRIBING  Tp  FLOW  IN  VORTEX  TUBES 


The  equations  of  motion  in  cylindrical  coordinates  for  a  three-dimensional 
laminar  vortex  flow  are 


a  (  >'  u  r) 

d  r 


^  a  (  r  u 

a  z 


0  Continuity 


(32) 


Ur  4-  U2 

a  r  a  z 


_L  i£  +  J,  -i-  fl  i.(rur  ) 

P  d  r  ar  [  r  ar 


+ 


Radial  Momentum 


(33) 


U  r  ^(f  U»)  du4> 

r  ^  r  ^  ^  z 


_a_  r_i_  ^ 

*'  ar  [r  ar 


a^u» 
a  z  ^ 


Tangential  Mcmentum 


(3M 


Ur 


aur 

ar 


+ 


Uz 


a  uz 

az 


_L  ^ 

az 


+  u 


-f  (' 

r  ar  ' 


Axial  Momentum 


(35) 


Although  the  flow  in  vortex  tubes  is  necessarily  viscous  everywhere,  the 
flow  in  the  main  vortex  (region  of  primary  flow  -  see  Fig.  l)  may  be  considered 
separately  from  the  thin  viscous  boundary  layer  in  the  end  walls  (region  of 
secondary  flow)  for  the  purpose  of  deriving  the  governing  equations  of  motion. 
This  assumption  is  justified  since  the  primary  flow  has  well-defined  and  neeirly 
constant  properties  in  the  axial  direction  and  any  strong  axial  gradients  are,  by 
definition,  confined  to  the  secondary  flow  region. 


Primary  Flow  Region 

The  important  properties  and  the  governing  equations  of  motion  for  the  pri¬ 
mary  flow  can  be  derived  from  Eqs.  (32)  through  (35)  by  an  order  of  magnitude 
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analysis.  Assume,  as  is  the  case  for  most  vortex  flows,  that  the  radial  velocity 
is  much  smaller  than  the  tangential  velocity,  i.e.. 


0  {€) 


€  <<  I 


(36) 


Then,  since  the  flow  is  injected  into  the  vortex  tube  with  constant  mass  flow 
per  unit  length,  the  order  of  magnitude  of  the  derivatives  in  the  z-direction  are 


dz 


and  the  order  of  magnitude  of  the  derivatives  in  the  r -direct ion  are 


au^ 

ar 


(37) 


(38) 


With  the  use  of  Eq.  (36),  the  zeroth  order  eqiiatlons  of  motion  for  the  primary 
flow  are 


a(rOr) 
a  r 


a  (r  Dz) 

■  "i  —  =  0  Continuity 


(39) 


P  dr  ^  ~ 


Radial  Momentum  (^O) 


+ 


1/ 


a^u» 

"az^ 


Tangential  Momentum 


(hi) 


az 


=  0 


Axial  Momentum  (^2) 
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Equations  (4o)  and  (42)  state  that  the  static  pressure  and  tangential  velocity  are 
a  function  of  r  only  and  not  a  function  of  z  .  Consequently,  the  radial  velocity 
is  GU.SO  a  function  only  of  r  .  This  result  has  been  verified  experimentally  in 
Ref.  10. 

To  satisfy  Eqs.  (39)  through  (42)  the  solution  for  the  primary  flow  must  have 
the  following  form: 


(r) 


(43) 


Ur  =  Ur  (  r  ) 


(44) 


P  =  P  (r) 


(45) 


Uj  =  f,(r  )  z  +  ) 


(46) 


This  class  of  solutions  with  an  important  generalization  for  a  finite  axial  pressure 
gradient  has  been  treated  extensively  in  Ref.  6.  By  introducing  Eqs.  (43)  through 
(46)  into  the  conservation  Eqs.  (39)  through  (42)  the  governing  equations  reduce 
to 


a(rur) 

dr 


-  r  f,  (r) 


1  ^ 

p  dr  ~  r 
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^  _^(ru,^)  .  j,  A 

r  dr  'dr 


j__d_(ru^) 
r  dr 


(^9) 


Laminar  Secondary  Flow  Region,  Including  Primary  Flow  Interaction 

The  secondary  flow,  illustrated  in  Figs .  1  and  2,  is  associated  with  the 
flow  in  the  boxandary  layers  at  the  end  walls  of  the  vortex  tube. 

Since  the  conditions  at  the  outer  edge  of  these  boundary  layers  are  determined 
by  the  primary  flow,  which  for  the  general  case  depends  upon  the  amount  of  secondary 
flow,  it  is  obvious  that  the  two  regions  cannot  be  analyzed  independently.  Conse¬ 
quently,  the  derivation  of  the  governing  equations  for  the  secondary  flow  included 
the  interaction  with  the  primary  flow. 

The  equations  governing  the  secondary  flow  were  derived  from  Eqs.  (32) 
through  (35)  by  integrating  over  Z  between  the  limits  0  <  Z  <  Z,  .  Two  basic 
assumptions  about  the  order  of  magnitude  of  terms  in  the  differential  equations 
were  made. 

=  0(0  (50) 


=  0  (€)  (51) 

I 


From  these  two  basic  assun^xtions  it  is  necessary  to  derive  several  auxiliary 
relations  to  simplify  the  equations.  Let  0^^  be  a  typical  tangential  velocity 
in  the  boundary  layer  and  Ur  a  typical  radial  velocity  in  the  boxindary  with 
unspecified  order  of  magnitude.  Then  in  the  boundary  layer  the  order  of  magnitude 
of  the  derivatives  are 


du<^ 

az 


(52) 
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d}^ 

dr 


0(%) 


t  =  °(^) 


(53) 


By  performing  the  order  of  magnitude  analysis  on  Eqs.(32)  through  (35)^  the 
following  resiiLts  are  obtained: 


!i'i 

Ur  8/ 


0(1  ) 


(54) 


(55) 


dz 


0(  c^) 


(^7^) 


(56) 


The  last  equation  is  significant  because  it  proves  that  the  axial  pressure 
gradient  in  the  end-wall  boundary  layer  may  be  neglected  in  a  zeroth  order  theory. 

The  order  of  magnitude  of  Ur  must  now  be  determined.  From  Eqs.  (33)^  (56) 
and  (4o),  it  can  be  shown  that. 


_LjLP 

P  d  ^ 


o  ■  0  (  -^r  ) 


(57) 


Then  with  the  use  of  Eq.  (55)  one  obtains  the  important  result  that 


=  0  (I) 


(58) 


Vi 


Finally,  by  substituting  Eq.  (58)  into  Eq.  (55)  one  obtains  the  result  that 


s_ 

r. 


O  ( - ! - \  =  0 


(€) 


(59) 
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The  order  of  magnitude  of  the  ratio  of  secondary  to  toteLL  flow  is  given  by  the 
parameter 


/3l 


27rr,  />  Ur  S 
27r  r,/)  Dr,  Z, 


2  Jj_ 

yRef  ilri  2 1 


(60) 


With  these  results,  if  terms  of  first  order  or  higher  in  €  are  neglected 
and  if  Eqs.  (47)  through  (49)  are  used  to  describe  the  primary  flow,  then  the 
governing  equations  reduce  to 


rur  Z, 


+  2  f  T  Uf  dz 

"o 


-n  Ur,  Z, 


Continuity 


(61) 


i: 


+  Uz 


d  Uf 

a  z 


r 


.  j.  ap 

P  dr 


-  V 


a  z 


dz 


Radial  Momentum 


(62) 


r*  r iii  4.  u,  a  _  y  a^uA  "I 

/  L  r  ar  ^  az  -  -aitj 


dz  =  0 


Tangential  Momentum 


(63) 


ap 

az 


=  0 


Axial  Momentum 


(64) 


For  convenience  it  is  also  determined  that 


(  r  Uj) 


z  : 


-/ 


dz 

ar 


(65) 


By  carrying  out  the  indicated  Integration  in  Eqs.  (62)  suid  (63)  two  integral 
boundary  layer  eqxxations  are  obtained. 
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o  o 

^  j (rur)(ru^)o'z  - 


ror^  d2  +  {  u/-  u/)  dz 


^  ^  dz  'z=o 


The  following  integral  quantities  are  defined 

=  27r/3j^  ru^  ru^c^z 


Moment  of  Moment um 


%  =  27rp  Tru^fl^z 


Secondary  Mass  Flow 


=  ZirpvT 


t^) 


/Z  =  0 

Torque  Per  Unit  Radius 


l  2  , 

=  ZTTpj  ru,dz 


Padial  Momentum 


r”  z  z 

P  =  ZTTpj  (u^  -  \i^)dz 


Pressure  Force  Per  Unit  Radius 


rr^  =  Zwpvr 


(^) 

\  dz  /z  =  0 


Shear  Force  Per  Unit  Radius 
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In  terms  of  these  newly  defined  quantities  the  Integral  equations  (57)  and 
(58)  become 


d  Mr 
d  r 


P 


-  rzr 

Radial  Momentum  Equation 


(74) 


dr 


Tangential  Momentum  Equation 


(75) 


Equation  (74)  states  that  the  rate  of  change  of  radial  momentum  is  equal  to 
the  difference  between  the  pressure  force  eind  the  radial  shear  force.  Equation 
(75)  states  that  the  rate  of  change  of  moment  of  momentum  is  equal  to  the  differ¬ 
ence  between  the  moment  of  momentum  brought  in  from  the  primary  flow  and  the 
torque  exerted  by  the  end  walls.  These  equations  are  equivalent  to  the  generEd.lzed 
three-dimensional  boundary  layer  equations  derived  in  Ref.  12. 

The  two  integral  boundary  layer  equations,  Eqs.  (74)  and  (75)^  contain  six 
unknowns  defined  by  Eqs.  (68)  through  (73)*  Therefore,  the  problem  now  becomes 
one  of  writing  four  additional  relations  among  the  six  unknowns.  Unfortunately, 
experimental  data  or  similarity  solutions  are  not  available  euid  therefore  it  is 
not  possible  at  this  time  to  obtain  the  necessary  relations  in  a  completely  satis¬ 
factory  manner.  In  view  of  these  difficiilties,  the  following  one-peurameter 
families  of  velocity  profiles  are  assumed  for  the  tangential  and  radial  velocity 
in  the  boundary  layer. 


F,  (-|-) 


(76) 


u 


r 


9i  R  3 


(77) 


where  Vp  is  defined  as  the  average  velocity  in  the  boundary  layer  by  the 
relation 
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(78) 


Finally  it  is  noted  that  for  laminar  flow,  one  boundary  layer  equation  can  be 
eliminated,  if  it  is  assumed  that  the  average  velocity  is  related  to  the  Karmen- 
Pohlhausen  parameter  (see  Ref.  12). 


U 


VT 


(79) 


By  substituting  Eqs.  (76)  and  (77)  into  Eqs.  (68),  (69),  and  (70),  the  following 
relations  are  obtained. 


M 


* 


(80) 


w. 


SL 


=  s)sjr  F^i/2 


(81) 


\  =  2-irpvT 


F. '  (0) 


(82) 


Substituting  these  three  equations  into  the  secondary  flow  tangential  momentiam 
equation,  Eq.  (66),  and  nondimensionalizing  the  term  yields 

^  [(RU^f  F3F,  dl  -  RU^  ^  [(RU^)"  ^^]J^F^dI-  |-RU^F;(0)  =0  (83) 

Similarly,  after  nondimensionalizing  the  radial  momentvun  and  tangential  momentvim 
equations  (Eqs.  (A8)  and  (49)  )  for  the  primary  flow  and  the  continuity  equation 
(Eq.  (61)  )  and  using  an  integrating  factor  for  Eq.  (83)  the  following  set  of 
eqmtions  is  obtained. 
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-■w  *  (w)'  /' 


A  =  {--A  i  ‘X} 


In  addition,  the  equation  for  the  static  pressure  coefficient  may  also  be  written 


P-  P, 


where 


F/(0) 


7 


Qi  F,F, 


i  (90) 

The  parameter  P^_  is  called  the  secondary  flow  parameter  and  is  a  measure  of 
the  fraction  of  the  total  flow  which  passes  through  the  boundary  layer.  When 

=  0(l)  a  significant  fraction  of  the  total  flow  moves  radially  inward  through 
the  end -wall  boundary  layer. 
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The  velocity  distribution  functions  F|  and  F2  were  chosen  as  simple 
polynomials  by  the  procedure  presented  in  Appendix  B  and  satisfy  the  following 
conditions : 


Z  =  0 


(91) 


2=1  F, 


F2  =  0 


f;  =  0  Fz  =  0 

F,”  =0  Fz  =  0 


(92) 


The  axial  velocitie 


i9  CUIU 


(65).  The  stream  function  is 


^  -- 


iU-f  —  — — 

R  U.  2  J  2 


(93) 


euid  the  equation  for  the  streamlines  is 


RU„ 


T 


(9^) 


The  coefficient,  ,  used  in  the  strong-vortex  solution  is  given  by 


(95) 
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Viscous  Core  Region 

In  the  viscous  core  region  (see  Fig.  l)  the  end -wall  boundary  layer  does  not 
exist  and  the  equations  describing  the  flow  are  Identical  with  those  for  the  pri¬ 
mary  flow  (Eqs.  (47),  (48),  and  (49)  ).  The  general  form  of  the  solutions  is 
given  by  Ey.  (46)  with 

f|(r)  =  -^  i  f2(r)  =  -b  0  <  r  <  (96) 

When  this  solution  is  inserted  into  the  continuity  equation  (Eq.  (47)  )  and  the 
result  matched  with  the  solution  for  the  primary  flow  the  following  relationship 
is  obtained: 


RUr 


(RUr), 


(97) 


Substitution  of  Eq.  (97)  into  the  tangential  momentum  equation  (Eq.  (84)  ) 
yields 


\ 

r  /  RCfUR,  R*  ^  1 

l«*P  1  R.  2  i  ■  'J 

[e.P  C"-';''  "•)-'] 

0  <  R  <  R*  (98) 


Thus,  the  solution  for  the  viscous  core  is  given  by  Eqs.  (97)  and  (98)^  provided 
values  of  X  ^  and  Ur,^  are  known  from  the  solution  for  the  primary  and  secondary 
flow. 


Turbulent  Secondary  Flow  Region 

When  the  primary  flow  vortex  is  sufficiently  strong,  the  suigular  momentim 
outside  the  boundary  layer  is  constant  ( strong -vortex  case)  and  the  secondary 
flow  may  be  analyzed  separately  from  the  primary  flow  and  the  viscous  core.  For 
this  case,  the  value  of  X  is  1  and  a  considerable  simplification  of  the  conserva¬ 
tion  equations  makes  it  possible  to  also  treat  the  case  of  turbulent  flow  in  the 
end-wall  boundary  layer. 
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The  turbiilent  velocity  distribution  is  chosen  to  be 


(99) 


Ur 


-  -  U^,  Vr  F2{-|-) 


(100) 


Equations  (99)  and  (lOO)  may  be  substituted  into  Eqs.  (68),  (69),  (70)  and 
(72)  to  obtain 


4> 


2  “  “ 

2  rr  p  r  Vr  S 


i: 


(101) 


^  =  2^^r,u^  V^Sj['Fcr(|)  =  K^ZTrr  ^  (^02) 


Mr  =  27r/,r,u^’'^V^^Sjr'F2V(|-)  (103) 

P  =  ^TT^JU^^SjTd  -  F|^)  (104) 


Finally,  in  the  absence  of  better  information  on  turbulent  end-wall  boundary  layer 
profiles,  the  shear  stress  equations  are  assumed  to  have  the  seime  form  as  for  a 
flat  plate - 

.0225 

^  Su^  (105) 
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.0225  ^ 

T, 


(lOC) 


The  velocity  Ur^  is  the  maximum  radial  velocity  in  the  boundary  layer  and 
is  the  axial  distance  at  which  it  occurs.  Substitution  of  Eq.  (105)  into 
Eq.  (70)  and  Eq.  (106)  into  Eq.  (73)  yields 


2v p 


.0225  0/ 

-  -  1/4 

(^) 


(107) 


''T, 


(108) 


Then,  substitution  of  Eqs.  (lOl)  through  (l04),  (107),  and  (108)  into  the  radial 
and  tangential  momentum  equations  (Eqs.  (7^)  and  (75)  )  give 


k3  ^(RVr^Ct)  = 


RVc 


(Ct  '^r) 


1/4 


(109) 


where 


ks  /r  (R  Vr  Ct) 
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E-,uations  (109)  and  (llO)  with  auxiliary  Eqs.  (ill)  through  (II9)  are  the 
required  equations  for  solution  of  the  secondary  flow  in  the  turbulent  end-wall 
boundary  layer . 


CONFIDENTIAL 


R -2494-1 


CONFIDENTIAL 


The  coefficient  KjCR)  used  in  the  solution  of  the  strong-vortex  case  i£ 
given  by 


K 


T 


(120; 


Boundary  Conditions  and  Method  of  Solution 

If  Eq.  (86)  is  substituted  into  Eq.  (85)  and  Eq.  (85)  is  substituted  into 
Eq.  (84),  a  single  second-order  differential  equation  in  X  results.  The  solu¬ 
tion  is  obtained  by  numerical  integration  with  the  boundary  conditions. 


R  =  1.0  X  =  1.0  Ct  =  0 


(121) 


R  =  R,^ 


C/\  . 


X*Re,  Ur, 


exp 


exp 


^  Rer  Ur, 

^  Re,  Ur,  I 


(122) 


The  numerical  solution  for  the  strong-vortex  case  was  obtained  by  setting 
;3l  =  0,  Rer  =  1000.  These  conditions  were  sufficient  to  maintain 


X  =  1.0  1,0  >  R  > 


(123} 


Furthermore  an  analytical  solution  for  the  strong-vortex  case  may  be  obtair.e 
by  substituting  Eq.  (123)  into  Eq.  (86).  The  result  is 


,  1/^  2/3  4/3  1/4 

(-f  a,)  R  (l-R  ) 


(I2t) 
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Equation  (124)  is  very  useful  in  comparing  the  effect  of  different 
['  velocity  profile  shapes  on  the  solution. 


1 

I 
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APPENDIX  B 


SELECTION  OF  THE  BOUNDARY  LAYER  VELOCITY  DISTRIBUTION  FUNCTIONS 


Laminar  Boiindary  Layer 

The  equations  describing  the  viscous  flow  adjacent  to  a  rotating  disc  have 
been  solved  exactly  (Ref.  lO).  As  a  means  of  selecting  suitable  velocity  profiles 
for  use  in  the  vortex  analysis,  the  radial  mass  flow  on  a  rotating  disc  was  eveiluated 
using  simple  polynomial  expressions  for  the  velocity  profiles.  Con^jarison  of  the 
approximate  mass  flow  computed  using  the  assumed  profiles  with  the  exact  mass  flow 
indicated  which  of  the  assumed  profiles  gave  the  most  accurate  mass  flow  rates. 

The  momentum  integral  equation,  Eq.  (66),  written  in  a  form  applicable  to 
the  rotating  disc  problem  becomes 


(I2p) 


since 


u^(0)  =  O 


(126) 


U^(8)  =  O 


(127) 


The  tangential  velocity  distribution  at  the  surface  of  the  rotating  disc  is  given 
by 


u,(0)  =  (ur 

V 


(12S) 


The  velocities  are  given  by 
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(f ) 


(129) 


(130) 


Then  substitution  of  Eqs .  (128),  (l29),  and  (l30)  into  Eq.  (l25)  yields, after 
some  mathematical  manipulation,  the  equation 


jr^  Fa  (0) 

^r\  ^ 


(131) 


where 


^  =  Vf  s 


(132) 


Since  77  is  a  constant  for  laimiuar  flow,  integration  of  Eq.  (131)  gives  the 
nondime ns ional  boundary  layer  thickness 


.  [i  1 


(133) 


The  induced  radial  mass  flow  at  r  is 


Q  =  Zvrju  dz  =  —  Zvr^^  vu>  f 
’'o  ^  ^  •^0 


(13A) 


Define  a  quantity,  B  ,  by  the  relation 


B  =  -  277’/'  F^dl 
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Hence  the  radial  mass  flow  becomes 

0  =  (136) 


The  velocity  distribution  functions,  F3  and  F4  ,  may  be  redefined  in  terms  of 
the  velocity  distribution  functions,  F|  and  F2  ,  employed  in  the  secondary  flow 
analysis  by  the  relations 

E3  =  I  -  P,  (137) 


F4  =  -*^2  (138) 

The  constant,  B  ,  which  depends  on  the  velocity  profiles,  is  given  in  Table  II  for 
comparison  with  the  exact  solution.  It  may  be  seen  from  the  results  presented  in 
this  table  that  velocity  profiles  and  yield  the  best  agreement  with  the 

exact  value  of  B  and,  consequently,  che  exact  secondary  mass  flow  rate.  Therefore, 
the  following  velocity  distribution  functions  have  been  employed  in  evaluating  the 
characteristics  of  a  vortex  tube  with  laminar  flow  in  the  end-wall  boundary  layer. 

F,  =  2Z  -  z*  (139) 


Fg  =  I2(z  -  2Z^  +  z’)  (l40) 


Turbulent  Boundary  Layer 

Since  a  simple  l/7-power  law  under  most  conditions  describes  flat-plate 
boundary  layers,  it  has  been  chosen  for  the  tangential  velocity  profile.  For 
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the  radial  velocity  profile,  the  l/7-power  law  was  also  chosen  with  the  addition 
of  a  polynomial  needed  to  satisfy  Eq.  (8l). 

F,  =  (f)  Ukl) 


These  distribution  functions  are  illustrated  in  Figs.  3  and  4. 
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LENGTH  OF  TUBE 


VELOCITIES  CONSIDERED  IN  ANALYSIS 
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END  WALL  OF  VORTEX  TUBE 


TANGENTIAL  VELOCITY  WITHIN  BOUNDARY  LAYER 
TANGENTIAL  VELOCITY  OUTSIDE  BOUNDARY  LAYER 
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VARIATION  OF  RADIAL  MASS  FLOW  WITH  REYNOLDS  NUMBER 

FOR  A  STRONG  VORTEX 


STRONG  VORTEX,  U^r  =  CONSTANT 
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FIG.  6 


VARIATION  WITH  RADIUS  OF  RADIAL  MASS  FLOW  WITHIN 
BOUNDARY  LAYER  ON  ONE  END  WALL  OF  VORTEX  TUBE 


STRONG  VORTEX,  Ci  ^  r  =  CONSTANT 
FOR  LAMINAR  FLOW,  WgL  =  27r^r|  EQ  (  6  ) 

FOR  TURBULENT  FLOW,  Wgj  =  2  7r;.i  f,  K  ^  Re  t°®  EQ  (  8  ) 
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DIMENSIONLESS  RADIUS, 
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EFFECT  OF  RADIUS  ON  ANGULAR  MOMENTUM  IN  PRIMARY  FLOW 
FOR  A  RADIAL  REYNOLDS  NUMBER  OF  10 
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EFFECT  OF  RADIUS  ON  STATIC  PRESSURE  COEFFICIENT 
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EFFECT  OF  RADIUS  ON  FRACTION  OF  TOTAL  FLOW  WITHIN 
BOUNDARY  LAYERS  ON  BOTH  END  WALLS  FOR  A  RADIAL 

REYNOLDS  NUMBER  OF  10 


'T'  AND  IN  FIGS.  3  AND  4 
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FLOW  STREAMLINES  FOR  RADIAL  REYNOLDS  NUMBER  OF  10  WITH  MODERATE 

SECONDARY  FLOW 
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tangential  velocity,  axis  of  vortex  tube 


FLOW  STREAMLINES  FOR  RADIAL  REYNOLDS  NUMBER  OF  10  WITH  LARGE 

SECONDARY  FLOW 
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TANGENTIAL  VELOCITY,  Ua  AXIS  OF  VORTEX  TUBE 
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FIG. 20 
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TANGENTIAL  VELOCITY,  AXIS  OF  VORTEX  TUBE 
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EFFECT  OF  RADIUS  ON  ANGULAR  MOMENTUM  IN  PRIMARY  FLOW 
FOR  A  RADIAL  REYNOLDS  NUMBER  OF  5 
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EFFECT  OF  RADIUS  ON  BOUNDARY  LAYER  THICKNESS  ON  END  WALL 
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EFFECT  OF  RADIUS  ON  ANGULAR  MOMENTUM  IN  PRIMARY  FLOW 
FOR  A  RADIAL  REYNOLDS  NUMBER  OF  20 
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